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Abstract 

We have extended the formalism for the two-neutrino double beta decay by includ- 
ing the weak-magnetism term, as well as other second-forbidden corrections. The weak 
magnetism diminishes the calculated half-lives in ~ 10%, independently of the nuclear 
structure. Numerical computations were performed within the pn-QRPA, for 76 Ge, 82 Se, 
100 Mo, 128 Te and 130 Te nuclei. No one of the second-forbidden corrections modifies signif- 
icantly the spectrum shapes. The total reduction in the calculated half lives varies from 
6% up to 32%, and strongly depend on the nuclear interaction in the particle-particle 
S — 1, T = channel. We conclude that the higher order effects in the weak Hamiltonian 
would hardly be observed in the two-neutrino double beta experiments. 



*Work partially supported by the CONICET from Argentina and the Universidad Na- 
cional de La Plata. 

tFellow of the CONICET from Argentina. 



PACS: 13.15+q;14.80;21.60Jz;23.40 



The sensitivity of the double-beta {(3(3) decay experiments is steadily and constantly 
increasing. To become aware of this it suffices to remember that, between the pioneer 
laboratory measurement on 82 Se |]J and the most recent on 76 Ge ||, the statistics has 
been improved by a factor of ~ 5000. 

The quantity that is used to discern experimentally between the ordinary Standard 
Model (SM) two-neutrino decays ((3(3 2v ) and the neutrinoless (3(3 events not included in 
the SM, both without {(3(3q v ) and with Majoron emissions (/3(3 M ), is the electron energy 
spectrum dY/de of the decay rate Y. It is usually given as a function of the sum (e) of 
the energies (ei and e 2 ) of the two emitted electrons. The (3(3 2v and (3(3 M decays exhibit 
continuous spectra in the interval 2 < e < Q, while the (3(3q v spectrum is just a line at 
the released energy Q = Ej — E F . 

In the evaluation of the (3(3 2v decay rate, the allowed (A) approximation is usually 
assumed to be valid, i.e., the Fermi (F) and Gamow- Teller (GT) operators are considered 
at the same level of approximation as in single- allowed (3 transitions. Besides, as the 
F operator is strongly suppressed in the (3(3 2v decay by the isospin symmetry, only the 
dominant contribution of the axial-vector current is considered in practice. Then, the 
(3(3 2v matrix element, between the initial state |0 7 ) in the (N, Z) nucleus to the final state 
|0 F ) in the (N — 2, Z + 2) nucleus (with energies E l and E F and spins and parities J n = + ) 
reads: 

wW, 2V (0+|^|1+) -(1+1^10+) 

N L l+ ~ L 0+ + Q/ 2 

Here the summation goes over the virtual states with spin and parity J 71 = 1 + . 

Except for the work of Williams and Haxton H, all higher order terms (HOT) in the 
weak Hamiltonian have been almost totally ignored in the past, simply because they are 
expected to be small. Yet, the question is not so simple. Firstly, from the comparison 
of the recent experimental results for the (3(3 half lives in 76 Ge 0: T 2v = 1.77 • 10 21 y, 
T M > 1.67 ■ 10 22 y and T 0u > 1.2 • 10 25 y, it can be stated that presently are being observed 
effects of the order of 10~ 4 at e ~ Q and of the order of 10 _1 at e ~ Q/2. Secondly, 
there are also several ongoing and planned experiments that are supposed to allow for 
measurements of still smaller effects. Thus, the HOT could, in principle, become relevant. 

In recent years we have examined the effects of the first-forbidden operators in the 
(3(3 2 u decay 0, |5|, both the non-unique (FFNU) and unique (FFU), which contribute via 
the virtual states J 71 = 0~, 1~ and J n = 2~, respectively. From these studies we have 
learned that: a) the FFNU transitions might increase the T 2v up to ~ 30%, yet they do 
not modify the A shape of the two-electron spectrum Q , b) the FFU transitions do alter 
the (3(3 2v spectrum shape but only at the level of 10~ 6 and mainly at low two-electron 
energy, where most backgrounds tend to dominate 0. Therefore, the effects of the first- 
forbidden transitions could hardly screen the detection of exotic (3/3 0u and (3(3 M decays, 
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which are the candidates for testing the physics beyond the SM. 

In the present work we want to complete the question of the competition between 
the standard and exotic /?/3-decays ||, by scrutinizing the effects of the nuclear matrix 
elements, which have been intensively studied in connection with deviations of the simple- 
beta spectra from the allowed shape || |7|, ||, |9]. (In particular, it should be reminded 
that the detection of the weak-magnetism (WM) term in the spectrum of /3 ± -decays 
of 12 N and 12 B to 12 C has provided a rather striking test of the CVC theory) It is 
customary to denominate these HOT as second-forbidden corrections (SFC), although the 
most significant among them comes from the WM term, which obeys the same selection 
rules as the GT operator. One also should keep in mind that the WM plays a very 
important role in both the (3j3 0l/ and the (3(3 M decays, through the so called recoil term 



Proceeding in the same way as in our previous works @, ^, [121 , we express the /3/3 2i . 
decay rate as: 



dY 



2v 



27ijj;\R 2u \ 2 5( 



ei + e 2 + uji + lu 2 



fc=i 



(2) 



where the symbol if represents both the summation on lepton spins, and the integration 
on neutrino momenta and electron directions. The transition amplitude reads [|], [5], |l"2]j : 



R2u 



2(2tt) 



£[1-P( ei e 2 )][l-P(^ 2 )] 



N 



(0+|g w (e 2 ^ 2 )|iv)(jv|g lv (e 1 z/ 1 )|0+) 
E N - Ej + e x + oJi 

(e^pi, s ei ) and v, 



(3) 



where the operator P exchanges the lepton quantum numbers e» 
(u>i, qi, s„J. All other notations have the usual meaning [[|. 

In the non-relativistic weak Hamiltonian we only retain the terms that are necessary 
for the evaluation of the SFC, and obtain 



H w (eiUi) 



Gg A 



[a + (3f + €i - qi ) X + q { (3f + 6i) Y + (3f + + 9i ) Z] ■ L(e^), (4) 



Here G 
aZ/2i? 



: (2.996 ± 0.002) x 10 12 is the Fermi coupling constant (in natural units), £ 
1.18Zy4 -1 / 3 is the Coulomb factor, and 



X 



9v 



(f w cr + r x p) 



3Mg A 

2 

r o- + 2 v / 8^(<T®F 2 (r))i 



r 

27 



(5) 



6M l 



2i(o- • p)r] 
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are the nuclear operators, with g v = 1 and f w = 4.7. The leptonic part is: 

1 „ , s +, , L cr ■ p 



L(eiz/j 



2e, ; 



F ( ei )x T (^) 1 - 



1 



<x(l - cr ■ q)x(-s^). 



(6) 



Keeping only the interference terms between cr with X, Y and Z, and the linear 
contributions in the lepton energies, from (§) and (f|), we obtain: Q 



R 



G 2 



2v 



12(2tt) 



/o(eie 2 )M^ + ^E^( e i e 2)-M 
2 i=i 



[l-P(^i^)]L(e 2 z/ 2 )-L(e 1 z/ 1 ), (7) 



where 

/o(eiea) 
and 



/i(eie 2 ) 



M 1 ' 



-MS 



-MS 



1 + 2 



ei + e 2 



/2(eie 2 ) 



ei + e 2 



6£-Q' J y ' Q 

(0+|<r|l+> - (1+|X|0+) 



/3(eie 2 ) 



2^(ee-Q)E 



TV ^l+-^0++ ( 2/ 2 

(0+M1+) • (l+|Y|0+> 



TV 



E 1+N -E 0j +Q/2 ' 
(0+|o-|l+) • (l+|Z|0+> 



(9) 



TV ~ + V/2 



It might be interesting to note that, while the FFNU transitions interfere with the GT 
operator at the level of the transition rates, the SFC do it already at the level of the 
transition amplitude. 

At the same order of approximation the differential transition rate reads: 



dY 



2/' 



4G 4 

157T 



i) 

2v- 



i=0 



where 



2v 



64tt 2 



(Q - ei - e 2 ) 5 J[ p k e k F (e k )de k . 



k=l 



Finally we derive the expressions for the spectrum shape 



(10) 



111 



(12) 



1 A factor of 3 has been omitted in the denominator of eq. (8) of ref. ||. All other formulas are, 
however, correct. 
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where 



Ti{e) = (Q- e) 5 f t (e) jT de lPl e 1 p 2 e 2 F (e 1 )F (e 2 ). (13) 
For the half-life we get 



T 2 M^+) = ln2[r 2 ,(0+^0+)]- 1 = [Mtj^QiM^ 



(14) 



with the kinematical factors 



As seen from (|12|) the spectrum shape mainly depends on the factors J-"i(e). They are 
displayed for 76 Ge in the upper panel of figure 1, as a function of the energy e. ^(e) an d 
•7-3 (e) exhibit the same energy dependence, while T\{e) shifts the A spectrum slightly to 
the right and ^(e) slightly to the left. 

The spectrum shapes without and with the SFC are compared in the lower panel of 
figure 1. It can be observed that they are quite similar. Analogous behavior was found 
for other experimentally interesting nuclei, such as 82 Se, 100 Mo, 128 Te, and 130 Te. 
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From (pi),® and flT0| ) one sees that the main effect of the WM consists in renormalizing 
the GT matrix element (H) as: 

M% - M® (l + ^) , (16) 

i.e., by a factor of ~ 1.05 for medium heavy nuclei, independently of the nuclear model 
employed. 

The matrix elements were evaluated within the pn-QRPA model, following the 
procedure adopted in our previous works 0, |12|, [L3| . We display them in table [[], together 



with the corresponding kinematical factors Gi- Besides the WM term, the velocity depen- 
dent matrix elements r x p and 2i(cr ■ p)r are also important, and particularly in the case 
of 100 Mo. The moment Ai^J is always relatively small. 



Table 1: Numerical results for the kinematical factors Gi (in units of y~ l ) and for the 
nuclear matrix elements Ai^l (i n natural units), evaluated within the pn-QRPA formalism 
with an effective axial charge g A = 1. 



Nucleus 


Go = G3 


Gi 


G2 


M$> 


M% 


M% 


M% 


76 Ge 


5.49 10~ 2U 


6.24 10- 2U 


2.33 10" 


-20 


0.050 


0.0044 


0.0002 


-0.0017 


82 Se 


1.83 1(T 18 


2.14 10~ 18 


0.83 10- 


-18 


0.060 


0.0041 


0.0003 


-0.0018 


100 Mo 


3.97 1(T 18 


4.54 10~ 18 


1.83 10- 


-18 


0.051 


0.0141 


0.0015 


0.0072 


128 Te 


3.54 10~ 22 


3.79 10" 22 


1.13 10- 


-22 


0.059 


0.0048 


0.0003 


-0.0016 


130 Te 


2.00 10- 18 


2.23 HT 18 


0.91 10" 


-18 


0.048 


0.0039 


0.0005 


-0.0014 



As it is well known, within the QRPA the GT moment JA^j strongly depends on 
the particle-particle coupling constant in the S — 1,T = channel, denoted by t in ref. 

(0) 
2v 



T3fl . This dependence is particularly pronounced in the physical region for £, where Ai^ 



goes to zero and the QRPA collapses. Ai^J an d Ai^J also strongly depend on t, but in 
a slightly different way than Ai^j- Thus the decay rates rely on two or three rapidly 
varying functions of t, which makes the calculated half-lives (T 2 ^, +,SFC ) to be even more 
sensitive on the value of t than within the A approximation (T£J). The numerical results 
are shown in table [2|. Contrarily to what happens in the case of the FFNU transitions, 
the SFC decrease the half-lives. The reduction ranges from 6% in 128 Te up to ~ 32% in 
100 Mo. 



2 From the theoretical point of view the nuclear matrix elements in 100 Mo are in same sense peculiar, 
because of the strong predominance of the [0g7/ 2 (n)0g9/ 2 (p); J x = 1 + ] configuration. 
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Table 2: Calculated half-lives within the allowed approximation (T^) and with second- 
forbidden corrections included (T 2 ^, +SFC ) in units of y. 



Nucleus 


rpA 

1 2u 


rjlA + SFC 

1 2v 


76 Ge 


7.3 10 21 


6.8 10 21 


S2 Se 


1.5 10 20 


1.4 10 20 


100 Mo 


9.7 10 19 


6.6 10 19 


128 Te 


8.1 10 23 


7.6 10 23 


130 Te 


2.2 10 20 


2.0 10 20 



This work completes our previous inquiries j|, |5| on the significance of the higher order 
corrections in the two- neutrino double beta decay. It can be concluded that: 

1) Both the first-forbidden transitions through the J 71 " = 0~, 1~ virtual states, and 
the second order corrections to the Gamow- Teller states, affect the transition rates in a 
significant way. 

2) The theoretical uncertainties within the QRPA, in the evaluation of the half-lives, 
are of the same order of magnitude (or even larger) than the experimental errors. 

3) The effect on the energy spectra of the higher order terms in the weak Hamiltonian is 
too small to shadow the possible exotic neutrinoless events in contemporary experiments. 
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Figure Caption 



Fig. 1. Kinematical factors jFj (upper panel) and the spectrum shape dr 2u /de (lower 
panel) for 76 Ge. All quantities are normalized to the maximum value of the allowed shape. 
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